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ABSTRACT 

A  study  and  pursuant  development  is  described  of  a  digital  computer 
program  for  the  computation  of  transfer  function  sensitivities  in  a 
symbolic  form  suitable  for  storage  and  subsequent  repetitive  numerical 
evaluation  over  a  range  of  frequencies  or  parameter  values.   The  algo- 
rithm is  implemented  in  conjunction  with  the  Network  Analysis  for  Systems 
Applications  Program  (NASAP)  version  developed  by  R.  S.  Schwartz  at 
Northeastern  University.   Several  example  problems,  as  well  as  sugges- 
tions for  possible  improvement  of  the  algorithm  and  its  extrapolation 
into  the  areas  of  optimization  and  automatic  design,  are  included. 
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I.   INTRODUCTION 

A.    STATEMENT  OF  THE  PROBLEM 

The  purpose  of  this  research  was  twofold:   (1)  to  investigate  the 
feasibility  of  providing  for  the  automatic  digital  machine  computation 
of  transfer  function  sensitivities  in  a  symbolic  form  suitable  for 
storage  and  subsequent  numerical  evaluation  over  a  range  of  frequencies 
and  parameter  values  or  for  use  by  a  minimization/optimization  procedure 
as  a  readily  available  source  of  gradient  data,  and  (2)  to  implement 
such  a  program,  if  feasible,  in  conjunction  with  the  Network  Analysis  for 
Systems  Applications  Program  (NASAP) . 

There  were  formidable  obstacles  encountered  in  the  course  of  the 
research  which  may  best  be  illustrated  by  a  preliminary  examination  of 
the  computer  core  storage  required  for  implementation  of  the  proposed 
procedure  if  it  were  meshed  with  an  existent  circuit  analysis  program 
such  as  NASAP.   The  NASAP  program  is  capable  of  manipulating  up  to  10,000 
independent  loop  products  (generated  through  use  of  topological  flowgraph 
techniques) ,  once  as  a  numerator  polynomial  in  s_  and  once  as  a  denominator 
polynomial;  that  is,  by  iterative  use  of  a  single  vector  variable,  LOOP, 
dimensioned  to  10,000.   Each  loop  product  may  represent  from  one  to 
twenty-three  constituent  elements,  and  the  analytical  expression  for  the 
partial  derivative  of  a  quotient  involves  four  separate  components  - 
numerator  polynomial,  denominator  polynomial,  partial  derivative  of  the 
numerator  and  partial  derivative  of  the  denominator.   The  implication,  then, 
is  a  storage  allocation  requirement  of  four  vector  variables  each  of 


dimension  230,000  words  (920K  bytes).   This,  of  course,  would  provide 
for  the  sensitivity  with  respect  to  a  single  parameter;  if  the  gradient 
is  desired  the  requirement  for  the  two  partial  derivative  variables 
would  be  multiplied  by  a  factor  of  23  again,  for  a  total  of  11,500,000 
words.   Consider  now  an  even  more  staggering  statistic.   Since  the  NASAP 
program,  like  most  programs  written  in  the  FORTRAN  language,  assigns 
only  one  alphanumeric  character  to  a  word  (four  bytes  in  the  IBM  360/67) , 
there  would  be  an  additional  factor  of  four;  three  for  the  element 
identifying  characters  and  one  for  a  code  to  indicate  the  algebraic 
formulation  required.   This  initial  examination,  then,  would  indicate  the 
need  for  an  impossibly  high  storage  allocation  of  46,000,000  words 
(184,000K  bytes!).   Some  degree  of  success  was  achieved  in  the  reduction 
of  this  figure  to  a  manageable  size,  but  at  the  expense  of  limitations 
on  the  original  objective.   The  extent  of  the  progress  and  associated 
restrictions  will  become  clear  in  later  sections. 

B.    DISCUSSION  OF  THE  NASAP  PROGRAM 

Prior  to  discussing  the  sensitivity  algorithm  in  any  detail,  it 
would  perhaps  benefit  the  reader  to  summarize  the  characteristics  and 
capabilities  of  NASAP  since  it  is  upon  this  foundation  that  the  proposed 
technique  is  firmly  based. 

The  revised  version  of  NASAP  is  a  digital  computer  program  for  the 
flowgraph  analysis  of  lumped,  linear  active  electrical  networks  designed 
to  provide  any  desired  circuit  transfer  function  in  either  symbolic  or 
numerical  form.   The  program  utilizes  the  "mixed"  method  to  generate 
circuit  equations,  i.  e. ,  operations  are  performed  on  a  tree-branch, 
link-branch  topology  so  that  Kirchoff 's  voltage  law  may  be  utilized  to 


express  each  link  voltage  in  terms  of  tree  voltages  and  Kirchoff 's 

current  law  may  be  used  to  express  each  tree  current  in  terms  of  link 

currents.   These  relationships,  in  association  with  volt-ampere  and  source 

dependency  relationships,  comprise,  a  formulation  defining  a  canonical  form 

called  a  primitive  flowgraph  [3]. 

The  user  of  NASAP  prepares  a  circuit  for  analysis  by:   (1)  converting 

to  a  linear  model  if  necessary,  (2)  labeling  each  circuit  element  by  type 

(R,L,C,V,I)  and  a  two  digit  number  which  uniquely  identifies  it,  and  (3) 

numbering  each  node  as  desired  except  for  the  ground  node  which  must  be 

designated  node  1.   The  user  then  inputs  this  connectivity  data  in  the 

following  format: 

ELEMENT     FROM       TO        ELEMENT     CONTROLLING  PARAMETER 
NAME      NODE      NODE        VALUE     (For  controlled  source) 

Following  this  are  coded  sets  of  data  describing  the  transfer  functions 
desired,  in  the  form  of  the  ratio  of  a  response  variable  (voltage  or 
current  associated  with  some  circuit  element)  to  an  excitation  variable 
(voltage  or  current),  and  the  forms  of  program  output  desired.   The  NASAP 
tree-link  hierarchy  [20]  and  the  order  in  which  the  elements  were  input 
to  the  program  determines  which  elements  are  placed  in  the  tree,  thus 
providing  the  user  with  some  flexibility  regarding  its  choice.-  Given  the 
topology,  the  machine  then  forms  the  fundamental  outset  matrix  and  from 
this  constructs  the  primitive  flowgraph. 

A  transfer  function  between  any  two  nodes  of  the  primitive  flowgraph 
may  be  obtained  by  adding  a  "closure"  [20]  and  solving  the  resulting 
"closed"  graph  by  application  of  Mason's  gain  formula  in  the  following 
manner.   If  the  transfer  function  between  node  j  and  node  k  of  a  flow- 
graph  is  equal  to  T.,  and  if  the  graph  is  then  closed  by  the  addition  of 


a  branch  from  k  to  j  whose  transmission  gain  is  equal  to  the  reciprocal 

of  T.,,  then  the  loop  gain  with  the  added  branch  is  equal  to  unity  and 

the  determinant  of  the  flowgraph  vanishes.   This  reduces  the  solution  of 

T   to  the  problem  of  summing  loop-set  transmittance  products,  a  procedure 
Jk 

highly  suitable  to  machine  solution.   Readers  desiring  more  detailed 
information  concerning  this  procedure  are  directed  to  Ref.  20. 

The  procedures  utilized  by  NASAP  for  labeling,  sorting  and  storing 
each  loop-set  consist  primarily  of: 

(1)  As  each  loop-set  is  determined  by  the  search  routine,  a  signed 
integer  label  is  generated  from  which  the  loop-set  may  be 
uniquely  regenerated. 

(2)  Loop-sets  with  the  same  label  but  opposite  sign  are  eliminated 
to  reduce  storage  requirements. 

(3)  Loop-sets  are  sorted  and  stored  in  order  of  decreasing  magnitude. 

(4)  The  completed  loop-set  array  represents  the  symbolic  solution  of 
either  the  numerator  or  denominator  polynomial  in  coded  form. 
From  successive  solutions  the  transfer  function  may  be  evaluated 
either  numerically  or  symbolically  by  decoding  and  summing  each 
loop. 

C.    A  COMPARISON  OF  TWO  FAVORABLE  TOPOLOGICAL  METHODS 

A  topological  method  other  than  the  flowgraph  procedure  just  described 
is  the  so-called  "k-tree"  approach  [6]  which  is  based  on  Maxwell's  topolo- 
gical  formulas  for  node  voltage  equations.   That  is,  transfer  functions  are 
derived  from  the  evaluation  of  determinants  formed  from  circuit  nodal 
equations ;however,  they  may  equivalently  be  formed  in  a  topological  fashion 
by  summing  all  tree  admittance  products  and  appropriate  2-tree,  3-tree,  ... 
k-tree  admittance  products. 
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In  both  methods,   k-tree  and  loop-set,  the  solution  of  passive 
circuits  is  obtained  as  a  sum  of  linear  products  of  the  network  parameters, 
thus  eliminating  the  requirement  for  operations  with  polynomials,  and 
facilitating  machine  evaluation  of  symbolic  solutions.   In  addition,  the 
two  methods  compare  favorably  with  respect  to  the  number  of  loop-trans- 
mittance  products  and  tree-admittance  products  for  a  given  circuit,  since 
the  maximum  number  of  trees  possible  for  a  given  topology  is  only  one 
greater  than  the  maximum  number  of  loop-sets  contained  in  the  corresponding 
primitive  flowgraph  [20].   Each  tree-admittance  product  is  unique;  however, 
in  the  flowgraph  method  several  loop-sets  are  frequently  found  to  have 
identical  transmittance  product  magnitudes,  although  only  one  will  remain 
after  cancellation.   Thus  each  remaining  product  is  unique.   A  primary 
advantage  of  both  procedures  is  this  feature  which  obviates  the  evaluation 
of  zero  terms,  as  they  would  be  if  the  system  determinant  were  evaluated 
by  matrix  algebra. 

Although  the  two  methods  appear  equal  from  the  above  discussion,  a 
very  important  difference  arises  in  the  solution  of  active  networks  since 
the  k-tree  approach  requires  the  graph  to  be  modified  and  divided  into 
separate  current  and  voltage  graphs  from  each  of  which  the  complete  k-trees 
must  determined,  whereas  in  the  loop-set  analysis  the  procedure  is  invariant 
for  passive  or  active  networks. 

Another  significant  disadvantage  of  the  k-tree  analysis  is  its 
dependence  upon  the  particular  transfer  function  to  be  evaluated,  i.e., 
impedance,  admittance,  gain,  etc.,  each  requiring  different  types  of  k-tree 
searches.   Thus,  the  number  of  k-tree  searches  will  vary  with  the  number  of 
transfer  functions  as  well  as  the  type  desired.   Therefore,  the  flowgraph 


method  appears  to  have  a  distinct  advantage  in  this  respect,  since  the 
procedure  remains  identical  regardless  of  the  number  or  types  of  transfer 
functions.   For  these  reasons  the  flowgraph  approach  was  chosen  for  imple- 
mentation in  NASAP. 

D.    SUMMARY 

It  was  the  objective  of  this  study  to  determine  the  feasibility  of 
utilizing  the  loop-set  symbolic  form  of  the  transfer  function  originated 
in  the  manner  described  in  (B)  above  as  a  basis  for  computing  and  storing 
sensitivity  information  in  similar  symbolic  form,  thereby  vastly  simplify- 
ing the  repetitive  computation  of  accurate  network,  sensitivities  and 
rendering  such  solutions  free  of  the  difficulties  involved  with  differen- 
tial, or  numerical  techniques. 


II.   DEVELOPMENTAL  PROCEDURE 

Utilizing  the  revised  version  of  NASAP  developed  for  the  PDP-10  and 
CDC-3300  computers  by  R.  S.  Schwartz  [20],  adaptation  to  the  IBM  360/67 
computer  at  the  Naval  Postgraduate  School,  Monterey  was  accomplished  by: 
(1)  altering  program  statements  necessary  to  reconcile  the  difference  in 
word  lengths  between  the  two  types  of  computers,  (2)  making  various 
syntax  changes  necessitated  by  technical  differences  in  the  implementa- 
tion of  FORTRAN  on  the  different  machines,  and  (3)  changing  the  modular 
form  and  adding  program  statements  necessary  to  provide  for  the  execution 
of  the  algorithm  as  intended  by  the  author.   Details  of  the  conversion  are 
apparent  from  comparison  of  the  program  listing  available  at  the  W.  R. 
Church  Computer  Center,  USNPGS ,  Monterey,  California  with  that  of  Ref .  20. 
The  listing  found  in  the  appendix  of  this  thesis  contains  only  the  NASI, 
NAS4,  and  SENS  subroutines.   The  first  two  are  included  since  extensive 
modifications  were  incorporated  in  the  course  of  implementing  the 
sensitivity  feature,  and  the  latter  because  it  is  entirely  new. 

As  previously  described,  the  NASAP  program  has  available,  in  symbolic 
form,  the  uncancelled  loop  products  created  from  a  topological  flowgraph 
analysis  of  a  given  circuit.   It  was  from  this  base  that  the  writer 
devised  a  program  to  compute  circuit  transfer  function  sensitivities  in 
an  algorithm  developed  as  follows. 

Since,  in  subroutine  NAS4  of  the  original  program,  there  implicitly 

exists  a  requested  transfer  function, 

T(s)  =  K  N(s)  in  which, 
D(s) 
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N(s)  =  bnsn  +  b,  sm  +  ...  +  b .  sP  +  ...  +  bT  n..s<l 
0      1  1  LPN 


D(s)  =  a„s   +  a,s   +  ...  +  a.s^  +  ...  +  a   s 
0       1  l  LPD 


q' 


where  n,  m,  p,  q,  n' ,m' ,  p',  q' ,  ...  are  integer  exponents  of  the  Laplace 
variable  s_,  not  sorted  into  any  ascending  or  descending  order  and  with 
repetitions  possible;  and  where  a.,  b.  represent  the  i   loop  products 
(algebraic  product/quotients  of  circuit  parameters  in  symbolic  form) ; 
and  where  LPD  is  the  number  of  uncancelled  loop  products  in  the  denomi- 
nator polynomial  of  T(s)  and  LPN  the  number  in  the  numerator.   Each  loop 
product  is  actually  determined  from  one  32-bit  computer  word,  wherein 
the  presence  (1)  or  absence  (0)  of  each  element  is  indicated  by  the 
corresponding  bit  value.   This,  of  course  limits  the  number  of  circuit 
parameters  which  can  be  accommodated  by  the  program  to  NBITS-1  =  31, 
although  in  the  current  version  this  upper  limit  is  set  to  23  due  to  other 
considerations . 

It  was  desired,  then,  to  compute  the  sensitivity  of  T(s)  in  symbolic 
form,  to  a  specified  circuit  parameter  P,  .   Since  sensitivity  is  essen- 
tially a  rate-of-change,  it  follows  that  the  transfer  function  sensitivity 

may  be  defined  as:  D(s).8N^s^  -  N(s) . 9D^S^ 

Sp  =  illsj_  ,  3(N(s)/D(s))    * k        ... 

P     3P         9  p 2 ^ 

k      k  k  D  (s) 

To  implement  this  fundamental  approach  utilizing  the  data  available 

in  NAS4,  the  following  simple  algorithm  was  developed,  (See  Flow-Chart, 

Fig.  1): 

(I)  Observe  the  j    factor  of  the  i   loop  product,  b    ,  in  the  numerator 


of  T. 


Is  b..  in  the  numerator  of  the  i   loop  product,  b.? 
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YES:     (a)  Is  this  factor  the  parameter  P,  ? 

YES:   (i)  Set  flag  to  indicate  that  b.  contains  P, 
b  1  k 


(b)  N.  .  =  b.  . 
ij     13 


NO:      (a)  Is  this  factor  P.  ? 
—  k 


YES 

:   (i)  Set  flag 

(ii)  Change  sign  of 

k'  1 

(iii)  Form: 

v   k'lj 

■  i/(b±j)2 

NO: 

Form  (*$.  \ 

=  1/b.. 
13 

3 

(b)  N.  .=  1/b.  . 
13     13 

(II)  Repeat  step  (I)  for  the  (j+1)    factor  of  the  ith  loop  product. 
Continue  until  all  factors  are  exhausted  (j>NFN(i),  where  NFN(i) 

represents  the  number  of  factors  in  the  i   loop  product  of  the 
numerator  polynomial) . 


ftl 


(III)    Does  (4S-    J     contain  Pfc?    (is    flag   set  to  1?) 


YES:      Continue 


■/3N   \ 


NO:        Set/il    ]       =     0 
Vi 


IV)  (A)    i   =   i  +  1 


(B)    Repeat    (I)    (II),    and    (III)    for   the    (i+l)th    loop  product. 
Continue  until  all   loops   are   exhausted    (i>LPN) . 
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(V)    Repeat  the  entire  procedure  beginning  with  step  (I)  for  the 

denominator  of  T(s). 
At  this  point  the  sensitivity  is  implicitly  available  in  the  desired  form 

(eqn.  1) ,  where: 

LPD  NFD(i)       EXPDS(i) 
D  =  i         n     D.  .  s 

i=l  j=l       « 

LPN  NFN(i)       EXPNS(i) 
N  =  z    n      N.  .s 
i=l  j=l      ij 

LPBNFD(i)  EXPDS(i) 


3P,     i=l  j=l     TF 
k        J      V  3  k  /  ij 


9N    LPN  NFN(i)  /    v      EXPNS(i) 
■22.     T.  H     (M  \     s 

3P 
k 


1=1   J  =  1     l^kiij 


■f-V* 

NFN(i)  =  Number  of  factors  in  the  i   loop  product  of  numerator 
polynomial 

NFD(i)  =  Number  of  factors  in  the  i   loop  product  of  denominator 
EXPNS(i)  =  Exponent  of  s_  in  the  i   loop  product  of  numerator 
EXPDS(i)  =  Exponent  of  s_  in  the  ith  loop  product  of  denominator 
One  feature  of  the  programming  which  became  necessary  in  order  to 
reduce  storage  requirements  for  the  symbolic  sensitivity  component  variables 
SNUM,  SDEN,  PARTLN,  and  PARTLD  was  a  packing  scheme  which  in  effect 
reduced  vector  sizes  by  a  factor  of  one-fourth.   In  the  NASAP  program,  as 
in  most  FORTRAN  programs,  each  alphanumeric  character  is  represented  by 
a  four-byte  word,  wasting  a  great  deal  of  core  storage  since  one  byte  is 
sufficient  for  one  character.   This  exemplifies  one  of  the  most  unattractive 
aspects  of  the  FORTRAN  language,  especially  when  dealing  with  character 
strings  or  literals  as  in  the  case  of  the  NASAP  symbolic  version. 

A  remedy  for  the  problem  which  was  utilized  in  this  program  is  a 
very  short  but  efficient  packing  subroutine  which  has  four  arguments: 


13 


(1)  a  four-byte  variable  to  be  "packed" 

(2)  an  integer  constant  specifying  the  byte  number  of  the  variable 
in  (1)  to  be  assigned  the  value  of  the  packing  variable 

(3)  packing  variable 

(4)  an  integer  constant  specifying  the  byte  number  of  the  variable 
in  (3)  which  is  to  be  assigned  to  the  variable  in  (1) 

Upon  calling  the  pack  subroutine,  dummy  variables,  locally  declared 
L0GICAL*1  (one  byte) ,  receive  data  passed  from  the  specified  byte  of  the 
packing  variable,  assign  this  to  the  specified  byte  of  the  variable  to  be 
packed,  and  then  return  this  new  "packed"  version  to  the  calling  program. 
In  this  manner  the  sensitivity  component  variables  were  packed  with  para- 
meter name  (e.g.  L01,  R04,  etc.)  in  bytes  one  through  three,  and  byte 
number  four  was  packed  with  a  code  number  to  indicate  the  form  in  which 
the  parameter  should  appear.   This  code  appears  below: 
0  —  Numerator  or  denominator  =  z,   partial  =  z 

3  —  Numerator  or  denominator  =  z,    partial  =  1 

4  —  Numerator  or  denominator  =  z,    partial  =  0 

6  —  Numerator  or  denominator  ■  1/z,  partial  =  1/z 

7  —  Numerator  or  denominator  =  1/z,  partial  =  -1/z^ 
9  —  Numerator  or  denominator  =  1/z,  partial  =  0 

With  this  scheme  the  program  was  rendered  capable  of  processing  net- 
works of  up  to  200  independent  loops  to  provide  for -sensitivity  with 
respect  to  a  single  parameter,  or  by  successive  solutions  to  provide  for 
any  or  all  sensitivities  desired.   With  this  loop-set  limit,  the  entire 
program,  including  NASAP  requires  180K  of  core.   A  guide  for  core  storage 
allocation  may  be  found  in  the  appendix  for  circuits  of  greater  complexity, 
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Subsequent  to  the  packing  operation,  data  is  transferred  to  a  new 
subroutine  SENS  for  manipulation  into  the  desired  form  and  for  numerical 
evaluation.   In  this  operation  it  was  found  impractical,  of  course,  to 
permit  each  product  element  to  occupy  a  separate  storage  location  in  the 
computer,  since  this  would  essentially  require  dimensions  on  the  order 
of  the  square  of  the  maximum  number  of  loop  products.   Therefore,  numerical 
evaluation  of  the  constituent  components  of  the  sensitivity  are  recursively 
computed  and  simultaneously  combined  at  this  point  in  the  program  to  form 
the  numerator  and  denominator  arrays,  PARTI  and  PART2.   Eliminating  all 
zero  vectors  and  combining  like  powers  of  s_  reduced  storage  requirements 
for  this  operation  essentially  to  twice  the  order  of  the  denominator 
polynomial. 
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1=1+1      LK1=L00P(I) 
LK2=LK1  KK=0 

KKK=0      FLAG=0      ISO 


N8=N8+1 


A 


Yes 


No 


•fejr 


N1=N8 
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IS=IS+IST(N1) 

NM=NAME(N1)/100 


No 


w 


KK=KK+1 
Transfer  NAM 
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Figure  la.   Flow  Diagram  of 
Sensitivity  Algorithm  in 
Subroutine  NAS4 
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No 
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Print 
Loop  Products 


Continue 
Processing 


Figure  lb .   Flow 
Diagram  of  Sensi- 
tivity Algorithm 
in  Subroutine  NAS4 
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1 
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No 

" 

PARTLD(I,J)= 
LNUM(J) 

FLAG=1 
PARTLD(I,J)=1 

\ 

1 
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LNUM(J) 
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Figure  lc.   Flow  Diagram 
of  Sensitivity  Algo- 
rithm in  Subroutine  NASA 
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1/LDEN(L) 


1 


FLAG=1 

PARTLD(I,J)= 

-1/LDEN(L)**2) 


<r 


DEN(I,J)= 
1/LDEN(L) 


Yes 


Print 
Loop  Products 


J=J-1 


-OH  4 


EXPDS(I)=IS 

SIGD(I)=JS 

NFD(I)=J 
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III.   USER'S  RULES 

Users  of  the  sensitivity  subroutine  must  be  familiar  with  the  simple 
NASAP  coding  rules  contained  in  Refs.  1  and  2  as  modified  below: 

(1)  Only  one  "OUTPUT"  and/or  "ROOTS"  transfer  function  request  should 
be  made  in  a  given  problem  since  the  sensitivity  subroutine  is 
designed  to  store  data  relevant  to  only  one  numerator  polynomial 
at  a  time.   Of  course,  the  user  may  wish  to  follow  one  set  of 
problem  cards  with  another. 

(2)  The  "SYMBOLIC  SOLUTION"  card  must  be  inserted  following  "EXECUTE". 

(3)  No  "MODIFY"  requests  are  permitted. 

(4)  The  card  following  "SYMBOLIC  SOLUTION"  must  contain  the  words 
(only  the  underlined  portions  are  essential  for  machine 
recognition) ; 

SENSITIVITY  PARAMETER=L02  (or  R01,  C09 ,  etc.) 
Any  number  of  blanks  or  identifying  characters  may  be  included 
so  long  as  the  underlined  portions  appear  somewhere  in  columns 
1-80  in  the  order  shown.   Notice  that  the  parameter  (with  respect 
to  which  the  sensitivity  is  desired)  is  indicated  by  the  consecu- 
tive characters,  (=) ,  (R,  L,  or  C) ,  (two  digit  element  number). 

(5)  If  it  is  desired  to  have  the  sensitivity  magnitude  and  phase 
evaluated,  the  next  card  in  sequence  must  contain  a  maximum 
frequency  in  columns  1-11  and  a  minimum  frequency  (not  zero)  in 
columns  21-31,  preferably  both  in  standard  FORTRAN  formal  Ell. 4, 
although  normal  floating-point  form  is  satisfactory  if_  all 
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significant  digits  (plus  sign)  remain  within  the  specified 
11  columns.   Specifying  these  frequencies  limits  will  result  in  a 
print-out  of  the  sensitivity  magnitude  and  phase  at  200  evenly 
distributed  frequencies  in  the  interval. 

(Caution;   If  the  user  omits  this  option  and  desires  to  follows 
the  current  problem  deck  with  others,  a  blank  card  must  be 
inserted  in  place  of  the  frequency  interval  card) . 
An  interested  user  may  provide  for  the  plotting  of  the  generated 
magnitude  and  phase  data  by  changing  the  variables  KMAGN,  ANGLE,  and 
FREQ  into  array  variables  in  the  SENS  subroutine,  then  passing  them  to 
an  existing  library  plotting  subroutine.    This  will,  of  course,  consider- 
ably increase  the  core  storage  required  for  the  program. 
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IV.   INTERACTIVE  APPLICATIONS 

The  NASAP  program,  in  conjunction  with  the  sensitivity  subroutine, 
may  be  easily  adapted  for  use  under  a  time-sharing  system  such  as  CP/CMS 
with  only  a  few  simple  statement  changes  which  depend  upon  the  user's 
choice  of  file  structure.   One  choice  successfully  used  at  the  Naval 
Postgraduate  School  is  to  remove  the  four  initial  statements:   (1)  CALL 
NASAP,  (2)  STOP,  (3)  END,  and  (4)  SUBROUTINE  NASAP.   It  is  possible,  then, 
to  "offline  read"  the  program  into  CP  in  the  following  modular  form: 

FILE  PROGRAM  CONTENT 

RSS  Program  RSS  (Driving  program) 

NASI  NASI,  FUNCTION  ITRAN 

NAS2  NAS2 

NAS3  NAS3 

NASA  NASA,  SUBROUTINE  PACK 

NAS5  NAS5 

SENS  SENS,  BLOCK  DATA 

Under  CP/CMS,  the  entire  program  can  then  be  implemented  by  executing 
"RSS",  with  data  being  input  at  the  terminal,  or  if  the  three  "READ" 
statements  in  NASI  are  changed  appropriately,  data  may  be  read  from  a 
separate  file  so  defined. 
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V.   DISCUSSION  AND  RECOMMENDATIONS  FOR  FURTHER  IMPROVEMENT 

As  the  reader  may  have  gathered,  the  feasibility  of  computing  sensi- 
tivities of  a  transfer  function  to  all  circuit  parameters  and  storing 
them  in  symbolic  form  appears  doubtful  under  present  storage  limitations. 
What  follows,  however,  is  a  recommendation  for  further  reducing  storage 
requirements  somewhat  while  maintaining  the  symbolic  concept.   If  only 
a  single  numerical  sensitivity  were  desired,  it  would  of  course  be 
considerably  more  efficient  to  recursively  evaluate  the  sensitivity 
from  successive  loop-sets  with  their  appropriate  values  first  as  given, 
then  in  partial  derivative  form.   This  procedure  would  eliminate  most 
of  the  storage  problems.   The  singular  advantage  of  the  symbolic  approach, 
however,  is  that  evaluation  of  the  sensitivity  function  at  various  fre- 
quencies, of  for  various  parameters,  is  a  very  simple  matter  of  substi- 
tution, rather  than  the  use  of  time-consuming  iterative  or  numerical 
techniques  which  may  pose  stability  or  accuracy  difficulties. 

It  is  apparent  that  substantial  changes  are  indicated  in  the  proposed 
algorithm  (or  perhaps  only  the  techniques  for  implementing  the  algorithm 
which  will  permit  the  sensitivity  subroutine  to  cope  with  circuits  of 
greater  complexity  by  means  of  a  reduction  in  the  large  storage  require- 
ments.  However,  it  must  be  remembered  that  the  resulting  exponential 
increase  in  the  number  of  independent  loops  with  respect  to  a  concomi- 
tant increase  in  the  number  of  nodes  [20]  will  in  itself  limit  the 
complexity  of  topological/f lowgraph  methods,  particularly  in  maximal- 
planar  networks.   Indeed,  the  difficulties  of  dimensionality  plague  all 
circuit  analysis  programs. 
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One  suggested  method  of  reducing  storage  requirements  would  be 
to  extend  the  previously  mentioned  concept  of  utilizing  individual 
bits  of  the  computer  word  to  represent  constituent  elements  of  the 
network.   Since  in  this  case  it  is  necessary  to  know  which  of  six 
mathematical  formulations  is  appropriate  for  each  element  of  the 
numerator  and  denominator  polynomials,  instead  of  only  two  formula- 
tions as  in  the  case  of  transfer  function  loops,  it  would  appear 
that  a  scheme  requiring  only  three  variables  (2^  =  8  possibilities) ,  for 
a  numerator  and  denominator  each  dimensioned  to  the  maximum  number  of 
loops  (10,000  in  the  present  version)  would  suffice  to  allow  sensitivi- 
ties to  be  computed  and  stored  in  symbolic  form.   This  would  repre- 
sent a  substantial  reduction  from  the  460,000  that  are  presently  required 
for  a  10,000-loop  circuit  to  accomplish  the  same  task.   This  suggestion, 
in  effect,  offers  an  execution  time-vs-storage  tradeoff,  since  a  coding/ 
decoding  algorithm  must  be  developed  which  will  intrinsically  involve 
a  tremendous  number  of  manipulations.   In  addition  it  would  require  an 
almost  total  rewrite  of  the  FORTRAN  programming,  particularly  that 
portion  meshed  with  the  NASAP  NAS4  subroutine,  and  it  may  very  well 
lead  to  difficulties  not  apparent  without  further  investigation. 
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APPENDIX  A  SAMPLE  PROBLEMS 


(I)  Sample  problem  No.  1  -  Emitter  Follower 


(a)   Emitter  follower  linear  model: 


(b)   NASAP  coding: 


NASAP  PROBLEM  NO.  1 
♦EMITTER  FOLLOWER 
VI  1  2  1 

11  1  4  1 
Rl  2  3  100 
R2  3  it  1M 
R3  3  1  1.5K 
Rij  I*  1  100 
CI  3  k    6PF 
C2  3  1-200PF 

12  k    1  0.0UVC2 

OUTPUT 

VR4/VV1 

EXECUTE 

SYMBOLIC  SOLUTION 

SENSITIVITY  PARAMETER=R04 
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(c)      Sensitivity   output   for  emitter   follower: 


NUMERATOR    POLYNOMIAL    :    SENSITIVITY    OF    VR0U/VV01      TO      POU 
0. 69UUE    kO    TIMES    S**       0 
0.1042E    31    TIMES    S**    -1 
0.66G7E    17    TIMES    S**    -2 


DENOMINATOR  POLYNOMIAL    :       SENSITIVITY 

0.2778E    19  TIMES    S**       2 

0.3473E    31  TIMES    S**       1 

0.1035E    hi  TIMES    S**       0 

0.1111E    30  TIMES    S**    -1 

0 .  28U5E    16  TIMES    S**    -2 
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(II)      Sample  Problem  No.    2   -  Butterworth  Filter 


(a)      Original   linear   circuit 


(b)      NASAP   coding: 


NASAP    PROBLEM 

♦THIRD    ORDER    BUTTERWORTH 

Rl    k   2    1 

R2    3    1    1 

CI    2    1    1 

LI    2    3    2 

C2    3    1    1 

VI    1    U 

OUTPUT 

VR2/VV1 

EXECUTE 

SYMBOLIC    SOLUTION 

SENSITIVITY    PARAMETER=R02 
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(c)      Sensitivity   output   for   Sample  Problem  No.    2   -  Butterworth 
Filter: 


NUMERATOR    POLYNOMIAL    :    SENSITIVITY    OF    VR02/VV01      TO      R02 
0.5000E    00    TIMES    S**       0 
0.5000E    00    TIMES    S**    -1 
0.2500E    00    TIMES    S**    -2 


DENOMINATOR  POLYNOMIAL    :       SENSITIVITY 

0.1000E    01  TIMES    S**       h 

0.U000E    01  TIMES    S**       3 

0.8000E    01  TIMES    S**       2 

0.1000E    02  TIMES    S**       1 

0.8000E    01  TIMES    S**       0 

O.^OOOE    01  TIMES    S**    -1 

0.1000E    01  TIMES    S**    -2 
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(Ill)   Sample  Problem  No.  3  -  Transformer  circuit 


(a)   Original  circuit  with  ideal  transformer: 


,25h    IDEAL 
l:-2 


5f 


=-"-0 


?3 


(b)   Linear  model: 


© 


(c)      NASAP   coding: 


NASAP  PROBLEM 

♦INPUT 

IMPEDANCE 

11 

1  2 

1 

Rl 

2  k 

3 

CI 

3  1 

0.5 

LI 

2  3 

0.25 

12 

3  1 

-2  IR1 

V2 

k    3 

-2  VC1 

OUTPUT 

REQUEST 

Vll/I  II 

EXECUTE 

SYMBOLIC    SOLUTION 

SENSITIVITY    PARAMETER=L01 

.1592E    03  .1592E    00 
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(d)      Sensitivity   output   for   Sample  problem  No.    3   -  Transformer 
circuit: 


NUMERATOR    POLYNOMIAL    :    SENSITIVITY    OF    VI 01/ I  I  01      TO       L01 

0.1000E    01    TIMES    S**       1 
-0.3333E    00    TIMES    S**       0 

0.2778E-01    TIMES    S**    -1 


DENOMINATOR  POLYNOMIAL    :       SENSITIVITY 

0.1778E    01  TIMES    S**       2 

0.2667E    01  TIMES    S**       1 

O.IU^E    01  TIMES    S**       0 

0.3333E    00  TIMES    S**    -1 

0.2778E-01  TIMES    S**    -2 
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(.e)      Sensitivity  magnitude   and  phase   output   for   Sample  Problem 
No.    3  -  Transformer   circuit: 

EVALUATION    OF    SENSITIVITY 


FREQUENCY 

MAGNITUDE 

PHASEON 

RA[ 

0.1592E 

00 

0.4352E  00 

0.1763E 

00 

0.958UE 

00 

0.9268E-01 

'  -0.1267E 

01 

0.1758E 

01 

0.5082E-01 

-0.1405E 

01 

0.2557E 

01 

0.3U98E-01 

-0.1U57E 

01 

0.3356E 

01 

0.2666E-01 

-0.1U8UE 

01 

0.U155E 

01 

0. 215UE-01 

-0.1501E 

01 

0.4954E 

01 

0.1806E-01 

-0.1512E 

01 

■0.575UE 

01 

0.1556E-01 

-0.1520E 

01 

0.6553E 

01 

0.1366E-01 

-0.1526E 

01 

0.7352E 

01 

0.1218E-01 

-0.1531E 

01 

0.8151E 

01 

0.1098E-01 

-0.1535E 

01 

0.8950E 

01 

0.1000E-01 

-0.1538E 

01 

0.b2bOE 

02 

0.1U32E-02 

-0.1566E 

01 

0.6330E 

02 

O.HaUE-02 

-0.1566E 

01 

0.6U09E 

02 

0.1397E-02 

-0.1566E 

01 

0.6489E 

02 

0.1380E-02 

-0.1566E 

01 

0.6569E 

02 

0.1363E-02 

-0.1566E 

01 

0.66U9E 

02 

0.1346E-02 

-0.1566E 

01 

0.6729E 

02 

0.1330E-02 

-0.1566E 

01 

0.6809E 

02 

0.1315E-02 

-0.1567E 

01 

0.6889E 

02 

0.1300E-02 

-0.1567E 

01 

0.B969E 

02 

0.1285E-02 

-0.1567E 

01 

0.1480E 

03 

0.60U9E-03 

-0.1569E 

01 

0.1U88E 

03 

0.6016E-03 

-0.1569E 

01 

0.1^96E 

03 

0.598^E-03 

-0.1569E 

01 

0.1504E 

03 

0.5952E-03 

-0.1569E 

01 

0.1512E 

03 

0.5921E-03 

-0.1569E 

01 

0.1520E 

03 

0.5890E-03 

-0.1569E 

01 

0.1528E 

03 

0.5859E-03 

-0.1569E 

01 

0.1536E 

03 

0.5828E-03 

-0.1569E 

01 

0.151+UE 

03 

0.5798E-03 

-0.1569E 

01 

0.1552E 

03 

0.5768E-03 

-0.1569E 

01 

0.1560E 

03 

0.5739E-03 

-0.1569E 

01 

0.1568E 

03 

0.5709E-03 

-0.1569E 

01 

0.1576E 

03 

0.5680E-03 

-0.1569E 

01 

0.1584E 

03 

0.5652E-03 

-0.1569E 

01 

0.1592E 

03 

0.5623E-03 

-0.1569E 

01 

ANS) 
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(IV)   Sample  Problem  No.  A  -  I.F.  Transformer 


(a)   I.F.  Transformer  circuit: 


(b)   Linear  model: 


R,f  C 


* 6 


©      ©      © 


■^W>  #  »W^     <l     c     o 


C,±  $R, 


© 


© 


(c)      NASAP   coding: 


NASAP    PROBLEM 
*I.F.    TRANSFORMER 
II    1    2    1 
Rl    2    1    .  5M 
CI    2    1    245PF 

R2    2    3    10 
LI    3    it    .U9MH 
L2    U    1    .01MH 
L3    Ij    5    ,J*9MH 
R3    5    6    10 
C2    6    1    2U0PF 
C3    2    6    2    PF 
Rit    6    1    ,5M 
OUTPUT 
VRU/I  II 
EXECUTE 

SYMBOLIC    SOLUTION- 
SENSITIVITY    PARAMETER=C03 
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(d)      Sensitivity   output   for   Sample  Problem  No.    4   -   I.F. 
Transformer:     ■ 


NUMERATOR    POLYNOMIAL    :  SENSITIVITY    OF    VROU/ 1  I  01       TO      C03 

0.ti720E    18    TIMES    S**  5 

0.4786E    16    TIMES    S**  It 

0.1293E    Ik    TIMES    S**  3 

0. U0U9E    10    TIMES    S**  2 

0.U725E    06    TIMES    S**  1 

0.2309E    02    TIMES    S**  0 
0.U000E-03    TIMES    S**    -1 


DENOMINATOR  POLYNOMIAL    :       SENSITIVITY 

0. U760E    kl  TIMES  S**  6 

0.4827E    kO  TIMES  S**  5 

0.1304E    38  TIMES  S**  k 

0.4095E    34  TIMES  S**  3 

0.U805E    30  TIMES  S**  2 

0.2375E    26  TIMES  S**  1 

0.U266E    21  TIMES  S**  0 

0.U07UE    15  TIMES  S**  -1 

0.1000E    09  TIMES  S**  -2 
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(V)   Sample  Problem  No.  5  -  Wien-Bridge  circuit 


(a)   Original  circuit 


o — =»» 


i o 


(b)   Linear  model: 


(c)   NASAP  coding: 


•- o 


V2  "  KVR. 


*— — — — o 


NASAP  PROBLEf! 

*WI EN-BRIDGE  CI RCU I T 

1112  1 

Rl  2  1  1 

CI  2  1  1 

R2  2  3  1 

C2  3  k    1 

V2  1  k    1  VR1 

OUTPUT 

VV2/I II 

EXECUTE 

SYMBOLIC  SOLUTION 

SENSITIVITY  PARAMETER=R02 
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(d)      Sensitivity   output   for   Sample  Problem  No.    5   -  Wien-Bridge 
circuit: 


NUMERATOR    POLYNOMIAL     :     SENSITIVITY    OF    VV02/II01       TO       R02 
IDENTICALLY    ZERO 


DENOMINATOR    POLYNOMIAL  :       SENSITIVITY 

0.1000E    01    TIMES    S**  k 

O.fcOOOE    01    TIMES    S**  3 

0.6000E    01    TIMES    S**  2 

0.U000E    01    TIMES    S**  1 

0.1000E    01    TIMES    S**  0 
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(VI)      Sample  Problem  No.    6   -  High-Q  Filter 


(b)      NASAP      Coding: 


NASAP    PROBLEM 

♦HIGH    -    Q    FILTER 

VI    1    2    1 

Rl    2    l»    20K 

R2    5    3    10K 

R3    6    it    20K 

RU    6    7    100 

R5    7    5    100 

R6    5    1    2K 

R7    7    1    1M 

CI    2    3    .02UF 

C2    5    it    .04UF 

C3    6    3    .02UF 

11  1    7    150    \Rk 

12  1    5    150    IR5 
OUTPUT 
VR6/VV1 
EXECUTE 

SYMBOLIC    SOLUTION 
SENSITIVITY    PARAMETER=C02 
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(c)      Sensitivity   output   for   Sample  Problem  No.    6   -  High-Q   filter 


NUMERATOR    POLYNOMIAL    :    SENSITIVITY    OF    VR06/VV01       TO      C02 

-0.2684E    09    TIMES    S**  6 

0.3125E    12    TIMES    S**  5 

0.3125E    09    TIMES    S**  h 

0.1031E    06    TIMES    S**  3 

,0.1250E    02    TIMES    S**  2 

0.5000E-03    TIMES    S**  1 


DENOMINATOR  POLYNOMIAL  :       SENSITIVITY 

0.1650E  24  TIMES  S**  6 

0.2869E  21  TIMES  S**  5 

0.1887E  18  TIMES  S**  k 

0.56U3E  Ik  TIMES  S**  3 

0.6908E  10  TIMES  S**  2 

0.1575E  06  TIMES  S**  1 

0.1000E  01  TIMES  S**  0 
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APPENDIX  B  GUIDE  FOR  THE  ALLOCATION  OF  CORE  STORAGE 


NUMBER  OF  LOOP-SETS  POSSIBLE 


(VARIABLE  DIMENSION) 

PROGRAM  CORE  REQUIRED 

*  200 

180K 

300 

205K 

400 

230K 

500 

255K 

600 

280K 

700 

305K 

**  800 

330K 

900 

355K 

1000 

380K 

***  5000 

505K 

10000 

630K 

** 


*** 


Corresponds  approximately  to  a  maximal  planar  network 

of  5  nodes 

Corresponds  approximately  to  a  maximal  planar  network 

of  6  nodes 

Corresponds  approximately  to  a  maximal  planar  network 

of  7  nodes 
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PROGRAM  LISTING 


66 


65 

300 

350 
2 

360 


351 
352 


1351 
1352 


1353 
1355 


1356 


1360 
1361 


1362 
1363 


1354 
1357 

3 

67 

68 

301 
303 


SU3R0UT 
INTEGER 
DIMENS I 
COMMON/ 
COM^DN/ 
,6), ITS 
( 15t6) , 
ITASK6( 
COMMON 
) ,LL,MO 
), 

NAM(23t 
M,LCONC 
,M2 

DATA  IB 
FRMAX=0 
FRMIN=0 
READ(5t 
DO  1  J  = 
IF( INLS 
DO  2  1  = 
IF( INLS 
IF (INLS 
CONTINU 
GOTO  3 
G0T0(4, 
CONTINU 
GO  TO  6 
GO  TO  ( 
IPRIN=2 
IND=4 
WRITE  ( 
WRITE  { 
READ(5, 
IFdNLS 
IF ( INLS 
SENSFL= 
DO  1353 
IF( INLS 
CONTINU 
KP1=<+1 
DO  1354 
IFdNLS 
PARAMC 1 
DO  13S0 
IF( INLS 
CONTINU 
PARAM(2 
DO  1362 
IF{  INLS 
CONTINU 
PARAMt  3 
READ(5, 
GO  TO  6 
CONTINU 
SENSFL= 
GO  TO  6 
K  =  J 

GO  TO  ( 
WRITE  ( 
WRITE  ( 
NASA=4 
WRITE  ( 
GO  TO  6 
GO  TO  ( 
WRITE  ( 


INE  NAS 
EOFCKF 
ON  INLS 
SEN/PAR 
DATA/I^ 
K3 

ITSK4(1 
17,12), 
IND,NAS 
DNM( 23) 

3) , IT  AG 
( 10,40) 

/'  • 

.0 

.0 

101  ,ENO 

1,32 

T(J )-IB 

1,3 

T(J )-I0 

T(J+1 J- 

E 


,  SENSFL,PARAM 

T( 80) ,N3DE(21) 

AM (3 ) ,SENSFL,FRMAX,FRMIN 

SRK  100)  ,  I  NSR 2(24,6)  ,ITSK1(  15,4)  ,ITSK2(3 

0,5)  ,ITSK5(9,10)  , I  RAY (19, 2) ,  ISYMB(48), 
I0PRTN(3, 2) 

A, MTRAN13) ,MD, IJTAG(15,2) , IPRIN,Z(42 
,VALM0D(23,10)  , I  LOOP (23) , L I , LIC , M9 , ZZ( 23 

(5,8) ,NMBR(6), 

, NAME (2  3,2)  ,MVI, I  ST (40) , INC (38,42) ,N1,NA 

/,IE/1HE/ 

=998)     INLST 

)300, 1,300 

PRTN(  1,1) )2, 350,2 
I0PRTN(I,2) )2,360,2 


301,351 ) ,1 

E 

6 

67, 303, 352, 68, 68), NASA 


6,119) 

6,  116HNLST 

101, END  =  1357  )     INLST 

T(1)-I0PRTN(3,1) ) 13  57,1351,1357 

T(2)-IE  J1357, 1352, 1357 

1 

K  =  3  PO 
T (K)-  IS YMB( 46 ) ) 1353, 1355, 1353 
E 


80 


.IE)  GO  TO  1354 

(I  ) 

0 

EO.ISYMB(MM)  )  GO  TO  1361 


0 


I=KP1, 
T{  I  ).E0 
)=INLST 

MM=1,1 
T(I+1). 
E 
)=MM-1 

MM=1,1 

T(I+2).E0.ISYMB(MM) )  GO  TO  1363 
E 

) =MM-1 
100,END=60)  FRMAX,FRMIN 

0 

p 

6 
o 

67, 69, 68, 68, 69), NASA 

6,119) 

6,109) 


6,116)1 
6 

67,303, 
6,304) 


NLST 

305, 68, 68), NASA 


38 


305 


400 


7 
69 


10 


13 


11 
12 

14 

15 
16 
17 

18 


308 
19 


GOTO  68 
K  =  J+1 
Kl  =  80 

1/  p  _  /L 

DO  400J=1,23 
ILOOP( J)=0 
MODNMC J J=0 


L2=1,1D 
J,L2)=0 


( 6,119) 

69 


DO  400 

VALMOD( 

Ll=6 

L2  =  l 

M9=l 

ISGN=1 

LD  =  0 

IZ=0 

LI1=0 

IS  =  1 

WRITE 

GOTO 

NASA=2 

Kl  =  30 

K2=0 

WRITE  (6 

WRITE  (6 

DO  5  1=1,21 

Z( I)=0. 

ZCI+21  )=0. 

NODE (I )=0 
DO  5  J=l,42 
INC( I, J  )=0 

MTRAN(1)=0 
MTRAN(2)=0 
MTRAN13  )=0 
MC  =  0 

N  =  0 

IPRIN=0 
M=2 

M0VFL=23 
ISGN=+1 
LD  =  0 
Ll=l 
M=M+1 

IF(M-M0VFL-3) 
L2=l 
IZ  =  0 
GO  TO 
WRITE 
DO  47 


119) 
116)INLST 


7,96,96 


66 

(6,116)INLST 

J=K,K1 
ILST  =  ITRAN(  INLSTl  J)  ) 
IF(  INSRKILST)  )  12,98,10 
L3  =  INSR1(  ILST} 
ILST=L3 

IFUNSR2CL3,  LI  1)99,98,  11 
L3=INSR2(L3,L1 )+l 
GO    TO    13 
L3  =  l 

ILST=-INSR1 (ILST)-l 
GOTO (14, 15,16, 17,18,3  0  8) , LI 
L3=ITSK1(L3,L2  ) 
1=1 

GO    TO    19 
L3=ITSK2(L3,L2  ) 
GO    TO    19 
L3=ITSK3(L3,L2) 
GO    TO    19 
L3=ITSK4(L3,L2) 
1  =  2 

GO    TO    19 
L3=ITSK5(L3,L2) 
GOTO    19 

L3=ITASK6(L3,L2) 
GO    TO    (20,21,21,24,26,27,28,29,30,31,33, 


39 


20 
21 

22 

23 

24 

25 

26 

27 

28 

29 

30 
31 


32 

33 
34 

35 

32  2 

36 

37 


38 

39 

375 

40 

41 
42 
43 
44 

310 

311 


36,37,3 

310,311 

375)  ,L3 

MN=10*M 

INC1MN+ 

IF(MN-N 

N=MN 

N0DE(MN 

G0T0(46 

IF(  ILST 

MN=ILST 

GO  TO  4 

ISGN=-1 

GO  TO  4 

INC( I,M 

GO  TO  4 

INC(I ,M 

GO  TO  ( 

I NC  ( I  •  M 

GO  TO  4 

LD=LD+1 

IZ=10*I 

IF(  ILST 

L2=L2+( 

GOTO  47 

IF{ INC{ 

LD=LD+1 

/100 

LD=LD+3 

IFCL1-6 

L2  =  6 

GOTO  47 

Ll=4 

L2=l 

ISGN=(4 

IF(  {  (IN 

ISGN=( 4 

MC=MC+1 

MTRAN(M 

MD=MC/2 

ME=MC-2 

GOTO  46 

Ll-5 

M2=M-1 

MC=1 

MT=1 

GOTO  7 

L1=L1+1 

L2=l 

GOTO  47 

WRITE  ( 

NASA=5 

JC=MC-2 

I J  TAG (M 

GOTD  46 

JC=MC-2 

IJTAGC M 

GOTO  (4 

JC=MC-2 

IJTAGIM 

GOTO  46 

MT=2 

MTRAN(  2 

GOTO  7 

MT=3 

MTRAN(3 

GOTO  7 

IS=-2 

G0T0(45 

MODNMCM 

LI  =1 

L2=4 


8,39, 40 t 4  5, 46, 47, 51 ,66,97,98,99, 41,42 ,43, 44, 
,312,313,314,316,316,319,316,316,316,316,321, 

N+ILST 

1 fM)=4* (L2-3)**3/ ( (L2-3 )**2+l)**2 

)  23,23,22 

)=NCDE( MN)+1 
,45,39) ,L3 
)99,47,25 


)=100*ISGN*ILST 

6 

)=INC(I ,M)+I SGN*MN 

39,46), I 

)=INC(I  , M)+ISGN*(10*MN+ILST  ) 

6 

Z+ISGN*ILST 
)99,47, 32 
120+(L2-K2)*(2-(L2-K2)*6) )/36 

l,M)/lD0-4)     34,35,34 

5+( ILST -10)* (582-(ILST-10)*(40+( IL ST- 10) *( 22) ) ) 

-ILST 

) 322,315,99 


-ILSD/3 

C(l ,M)/100)-4)/3)46,98,46 

-ILST)/3 

T)=MD 

*MD+1 


6,376) 

*MD+1 
D,JC)=100*ISGN*ILST 

*MD+1 

D,JC)  =1  JTAGC-MD,  JC)  +  ISGN*MN 

5, 46), ME 

*MD+1 

D,JC)=I JTAGCMD, JC ) + I SGN*( 10*MN+ILST ) 


)=MTRAM(1) 


)=MTRAN( 1) 


,46),L2 
9)=100*ILST 


40 


GOTO   47 

312  MODNM( M9) =MODNM( M9 ) +1 0*MN  + I LST 
GOTO    46 

313  MODNMC M9)=M0DNM[ M9)+MN 
GOTO    45 

319  IS= IS+1 

IF( IS* LI +3) 99, 3 15, 3 24 

314  IF(LI-10)315,315,326 

315  V  =  FLOATUZ) 

VALMODl  M9,LI  )=V/10.0**  (LD) 

L2  =  12 

IF( IS)331,99,332 

331  IF{LI-3)47,334,324 

334  IF(L3-36)47,340,47 

340   IL00P(M9  )  =  (VAL'40D{M9,1  )-V  ALMD  D(  M9,  2  )  )  /V  ALMOD(  M9,  3  )-  1.  0 
10  001 
IF( IL00P(M9) ) 333, 32  8, 3  28 

332  IL00P(M9)=LI 

33  3   IF(LI1-IABS(IL00P( M9)  ) ) 335,47,47 

335  LI1=IABS(IL00P(M9) ) 
GOTO  47 

316  IZ=0 
LD=0 
ISGN=1 

IF(L3-35) 317,317,323 
323   LI=LI+1 
L2=7 

I  3=  I   3—  "\  (~\ 
G0T0(31,26,45,47),L3 

317  IF(IABS(IS)-1)99,338,324 

338  IS=1 
L3=L3-33 
GOTO(320,337),L3 

320  IS=-2 
L2=3 

33  7      M9=M9+1 

IF(M9-M)339,339, 326 

339  G0T0(47,311) ,L3 

321  IF(M0DNM(M9) / 100-4) 34,3  5,34 

45  L2=L2+1 

46  L2=L2+1 

47  CONTINUE 

IFtLl-6  )330,336,99 

336  LI=LI1 
IND=4 

IF(IABS( IS)-1)99,60,324 
330       IF( ( 10*Ll+L2)/ll-3)98,48,49 

48  INC(2,M)=0 

49  1=1 
IF(Ll-5)50,7,99 

50  Z(M)=FLOAT( IZ) 
Z(M)=ZIM)/10.D:fc* (LD) 
GOTO    6 

51  N1=N+1 
IF(MTRANd))     99,95,202 

20  2       IF( IMLST{J+8)-ISYMB(26) ) 204,52 ,2  04 
204       IF(INLST(J+8)-ISYMB(13 ) ) 5 5, 20 1,55 
201       DO    203    J=l,2 

DO    200    1=1,19 

K=IRAY( I, J ) 
203   IRAY( I , J)=ISYM3( K) 
203   WRITE  (6,116)(  IRAY( I ,J)  ,1  =  1, 19) 

52  IPRIN=1 

55  NA=0 

DO  58  1=1, N,l 

IF(NODE(  I)-l )  58,56,57 

56  WRITE  (6,112)  I 

57  NA=NA+1 
NODE(NA)=I 

58  CONTINUE 
IF(NASA-2)70,71,70 

70   NASA=4 


41 


GO    TO    66 

71  NASA=3 
IND=2 

IF(N-NA)     99,60,59 

59  WRITE     (6,111)     (NODE(I),     I=1,NA,1) 
WRITE     (6,110) 

GO    TO    60 
91  IND=7 

60  RETURN 

95  WRITE     (6,115) 
GO    TO    90 

96  WRITE     (6,114)     MOVFL 
90       NASA=4 

GOTO    66 

97  WRITE  (6,106) 
WRITE  (6,116)  INLST 

72  NASA=5 
GO  TO  49 

98  WRITE  (6,108)  J 
WRITE  (6,116)  INLST 
IF(l_l-6)72,66,  99 

99  WRITE  (6,107 ) 
GO  TO  90 

326  WRITE  (6,327) 
GOTO  66 

328  WRITE  (6,329) 

324  WRITE  (6,325) 
GOTO  66 

100  FORMAT( Ell. 3, 9X, Ell. 4) 

101  F0RMAT(80A1) 

106  F0RMAT(24H    FUNCTIONS    NOT    AVAILABLE) 

107  F0RMAT(27H    MACHINE    LANGUAGE    BREAKDOWN) 

108  F0RMAT(32H  INPUT    CODING    ERROR    IN    COLUMN, 13) 

109  F0RMAT(39H       USE    HEADING    LABELED  NASAP    PROBLEM) 

110  F0RMAT(/,1X,33HTHE    ANALYSIS    WILL    CONTINUE    ON    THE    ASSUM 

1  ,5HPTI0N, 

2  26H    THAT    THE     LIST    IS    CORRECT. ,/, IX ,13HH0WEVER,     SKIP 
3,9HPING    NODE, 

4         39H    NUMBERS     IS    NOT    A    RECOMMENDED    PRACTICE.//) 

111  FORMATt 1X,41HTHE    FOLLOWING    NODE    NUMBERS    ARE    MENTIONED 
1  ,6HIN    THE, 

212H    INPUT    LIST     , / , IX, I  2 , 19 ( 1H , , I  2 , IX ) ) 

112  FORMATt//, IX, 5HN0DE    ,I2,27H    HAS    ONLY    ONE    ELEMENT    CONNE 
111HCTED    TO    IT. ,//) 

114  FORMAT( 1X,12HN0    MORE    THAN,I3,18H    ELEMENTS    ALLOWED.) 

115  FORMATt IX, 24HTHERE     IS    NO    OUTPUT    LIST.) 
116       FORMAT (1X,80A1  ) 

117    FORMAT( IX, 1116) 
119       FORMAT  (1H1) 

304       FORMAT( //l X, 39 HMODI F Y    AND    SYMBOLIC    SOLUTION    MUST     FOLLO 
1,9HW    EXECUTE// ) 

325  F0RMAT(//1X,39HF0R    LOOPS    YOU    MUST    SPECIFY — (START, END, 
1,4HINC) ) 

327  F0RMAT(//1X,39HT00    MANY    ELEMENTS    OR    VALUES    IN    MODIFY    S 
1  ,8HTATEMENT) 

329  F0RMAT(//1X,39HSTART    GREATER    THAN    END    AND/OR    INC    TOO    L 
1 ,4HARGE) 

376       FORMAT    (/1X,27HTHE    TRANSFER    FUNCTION    INPUT/ 
+1X,33HVARIABLE     IS    A    DEPENDENT    QUANTITY.) 
998    STOP 
END 


FUNCTION    ITRANU  ) 

ITRAN=1+I/16777216 

IF(  ITRAN)1,2,3 

ITRAN=ITRAN+63 

RETURN 

ITRAN=ITRAN-31 

RETURN 

WRITE(6,4) 
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RETURN 

FORMAT( //1X,15HMACHINE    FAILURE) 

STOP 

END 


SUBROUTINE    NAS4 

INTEGER    SENSFL,FLAG,EXPNS,SIGN,WHERE,EXPDS,SIGD, 
1PARAM,PARTLN,PARTLD 

DIMENSION    J  1(5)  ,  Ml  5)  ,  LNUM  (  23,  3  ) ,  LDEN  (  23  ,  3  ) ,  ILGTH(23)  , 
1IFMTK  18)  ,IFMT2( 12) 

COMMON /SEN/ PARAMO  J  ,SENSFL,FRMAX,FRMIN 

COMMON    IND,NASA,MTRN1,MTRN2,MTRN3,MD, IJTAG(15,2),I 
1PRIN,Z(42)  ,LL,M0DNM(23)  , VALMOD( 23 , 10)  ,1  LOOP (23)  ,LI,LIC 
2 ,M9, ZZ(23) , NAM (2  3, 3 ), ITAG( 5,8) ,NMBR(6), 
3M,SNUM(3  0,5)  ,SDE N( 3  0 ) , I OD (30 )  , I  0 ( 5 ) , NUM ( 5 ) 
5,XG(5),IP(30,5) ,  ID,NUMZ,DMY( 23) 

4,  NAME  (23  ,2  )  t  M\/  1 ,  1ST  (  40  ) ,  L  IM  (  6  ) ,  I  TWO  (23)  ,LOOP(  10000) 
6, WHERE , PART LN(  2  00,2  3)  tEXPNS 
7(200),NFN( 200)  ,PARTLD<  200,23) , EXPOS (200)  , 
8SIGDO00)  ,NFD(200)  ,  LPN  ,  LP  D,  S  I  GN  (  200  ) 

DATA     ILSTH/1H1, 1H2, 1H3 , 1H4 , 1 H5 , 1H6 , 1H7 , 1H8 , 1H9 , 2H10 , 
12Hll,2H12,2H13,2H14,2H15,2H16,2H17,2H18,2H19,2H20t2H21 
?     P  H  ?  ?     P  H  ?  3  / 

DATA     IFMT1/4H( 1X,,4HI5,I,4H8,3X,4H,A1,,4H1H(,,1H1, 
14H(A1,  ,1H2,4HI 1),,4H3H    / , 2H    , ,  1H1, 4H( A  1 , , 4H2I 1 ) , 
24H,6H),4H    S    * , 4H* , I  3 , 1H ) / 

DATA    IFMT2/4HC  IX , , 4HI 5 , I ,4H8, 3X ,4H , Al , ,4H 1H(  ,  ,  1H1  , 
14H(A1, ,1H2,4HI 1) ,,4H6H)     ,4HS    **,4H,I3)/ 

DATA    IM/«-  '/,IB/»  '/,IA0C/'A1,    '/,lTVl/'2  '/ 

MDD=MTRN1+1 

IF ( I PRIN-1)1,1 ,23 

1  IERR=0 
IF(LIC)60, 2, 5 

2  DO  3  Nl=l, MVI 

3  ZZ(N1)=Z(N1) 

IF( I  LOOP (1 )-10  00)5, 16,60 
16    DO  19  Nl=l ,MTRN1 

DO  19  L=l,  2 
19    IJTAG(N1,L)  =  IJTAG(N1,L  J-1000* ( I JTAG( N 1, L ) / 1000) 

IF( IPRIN)60,13,23 
23    1  =  0 

LPN=0 

DO  1024  K2=1,MDD 

N7=LIM(K2) 

WHERE=0 

IF(K2-MDD) 50,51,60 

50  WRITE  (6,135)  (  ITAGC K2 , K ) , K= 1, 8 ) 
WRITE  (6,160) 

J  =  -l 

IF(NMBR(K2) )601,52,600 

600  JS=IM 
GOTO  602 

601  JS=IB 

602   WRITE  (6,162)  JS 
GOTO  52 

51  WRITE  (6,138) 
WRITE  (6,160) 
JS  =  IB 

WRITE  (6,162)  JS 

SIGD( 1)=JS 

NFD(1)=1 

LPD  =  1 

J=l 

52  N8=l 
IF(N7)60,55,55 

55  WRITE  (6,100) 
WRITE  (6,141) 
GO  TO  1024 

56  DO  24  N9=1,N7 


43 


44 
25 
26 
27 
30 

28 


31 

230 

32 

231 

33 

29 

499 
500 
502 
509 

510 


511 

519 

530 
531 

534 


532 
533 


54  0 

542 
554 
555 

799 


800 


FLAG=0 

KK  =  0 

KKK=0 

ISLOP=0 

JS  =  IB 

1  =  1+1 

LK1=L0C 

LK2=IAB 

IF(J*LK 

JS=IM 

NL=N8 

IF(LK2~ 

IF(N1-N 

N8=N8+1 

GOTO    25 

ISLOP=I 

NM=NAME 

IF( (NM+ 

KKK=KKK 

DO    230 

LDENKKK 

GO    TO    3 

KK=KK+1 

DO    231 

LNUMIKK 

LK2=LK2 

N1=N1+1 

IFtNl-M 

IF(SENS 

JJ  =  1 

IF(K2-M 

LPN=LPN 

IFtLPN. 

IFtKK.E 

JJJ  =  1 

DO    511 

IFCK.GT 

CALL    PA 

CONTINU 

IF(LNUM 

1ND.LNUM 
CALL    PA 
GO    TO    5 
FLAG=1 
CALL    PA 
JJ=JJ+1 
IF( JJ-K 
IFCKKK. 
LX  =  0 
JJ=JJ-1 
JJ=JJ+1 
LX=LX+1 
JJJ  =  1 
DO    532 
IFCK.GT 
CALL    PA 
CONTINU 
IF(LDEN 

1ND.LDEN 
CALL  PA 
GO  TO  5 
FLAG=1 
CALL  PA 
IF(LX-K 
JJ=JJ-1 
WHERE=1 
GO  TO  3 
LPD=LPD 
IF(LPD. 
IFCKK.E 
JJJ  =  1 
DO    812 


PCI) 

SCLK1) 
1)44,60 


ITW0(N1 
8  )60,3D 


,25 


) )27,28,28 
,29 


SLOP+IS 
(Nl,l)/ 
3)*(NM- 
+  1 

K  —  1      3 

K,K)=NAM(N1,K) 
3 


TCN1J 

4)*(NM-5) ) 32, 31, 32 


K  =  l,3 
,K) =NAM 
-ITWOM 

VI-D26 
FL-1)    3 


(N1,K) 
1) 

,499,60 
4, 500,60 


DD) 509,799,60 
+  1 

GT.200) 
0.0)    GD 


WRITE(6,164) 
TO    531 


K=l,3 

.1)    JJJ 

CK(PARTLN(LPN,JJ),K,LNUM( JJ,K) ,JJJ) 

E 

( JJ,1). 

(JJ,3). 

CKtPART 

30 


=  4 


EO.  PAR  AMID.  AND.  LNUM(  J  J,  2)  .  E  Q.  PAR  AM  (2)  .A 
E0.PARAM(3))     GO    TO    519 
LN(LPN,JJ) ,4,0,4) 


CK(PARTLN(LPN,JJ)  ,4,3,4) 


K-D510 
EQ.O)    3 


,531,60 
0    TO    554 


K  =  l,3 

.1)    JJJ=4 

CK(PARTLN(LPN, J  J ) , K , LDEN ( LX , K ) , JJJ ) 


(LX,1  ). 
( LX,3). 
CK(PART 
42 

CKtPART 
KK) 534, 


EQ.PARAM(1).AND.LDEN(LX,2) . EQ.PARAM  (  2 )  .  A 
E0.PARAM13) )     GO    TO    540 
LN(LPN,JJ) ,4,6,4) 


LN(LPN,JJ) ,4,7,4) 
555,63 


4 
+  1 

GT.200  ) 
0.0)    GO 

K  =  l,3 


WRITE(6, 164) 
TO    831 


44 


JJJ=4 


IF(K.GT.l)     JJJ 

CALL    PACKtPART 
812       CONTINUE 

IF(LNUM( JJf 1). 
IND.LNUMt JJ,3) . 

CALL    PACKtPART 

GO    TO    830 
819       FLAG=1 

CALL    PACKtPART 

830  JJ=JJ+1 

IF( JJ-KK-1 )800 

831  IF(KKK.EO.O)     3 
LX  =  0 

JJ=JJ^1 
834   JJ=JJ+1 
LX=LX+1 
DO  832  K=l ,3 
JJJ=1 
IFtK.GT  .1) 
CALL  PACKtPART 

832  CONTINUE 

833  IF(LDENtLX,l  ). 
IND.LDENt  LX,3). 

CALL    PACKtPART 

GO    TO    842 
840      FLAG=1 

CALL    PACKtPART 
842       IFtLX-KKK)834, 

854  JJ=JJ-1 

855  WHERE=2 

34   IF(KK)60,42,43 

42  KK=1 
LNUMtltl  )=IB 
LNUMtl»2)=IB 
LNUMd,  3)=1 
IFMT1(8)-IA0C 

43  IF(KKK)60,40,4 

40  IFMT2(6)=ILGTH 
WRITEt6,IFMT2) 

1L0P 
GO    TO    901 

41  IFMT1(6)=ILGTH 
IFMTlt 12)=ILGT 
WRITE16,  IFMT1  ) 

1LDEN1K,K1) ,K1= 
IF{ IFMT1C8)-IT 
232        IFMT1(8)=ITW 

901  IF(  WHERE-D24, 

902  EXP^S(LPN)=ISL 
SIGN(LPN)=JS 
NFNt  LPN)=JJ 
IF(FLAG)60,905 

904  EXPDS(LPD)=ISL 
SIGDtLPD)=JS 
NFDtLPD)=JJ 
IF(FLAG)60,90S 

905  DO  906  K=ltJJ 
ITEST=0 

CALL  PACKt ITES 
IF { ITEST-5)960 

960  CALL  PACKtPART 
GO  TO  906 

961  CALL  PACKtPART 

906  CONTINUE 
GO  TO  932 

908   DO  931  K=1,JJ 
ITEST=0 

CALL  PACKtITES 
IF { ITEST-5J962 

962  CALL  PACKtPART 
GO  TO  931 

963  CALL  PACKtPART 


=4 

LDtLPD,JJ),K,LNUM(JJ,K) , JJJ) 

EQ.PARAM(l) . AND.LNUMt J J, 2) . EQ. PAR AM (2 ) . A 
E0.PARAMt3))     GO    TO    819 
LDtLPDt JJ ) ,4tO,4 ) 


LD(LPD,JJ) ,4,3,4) 

,831,60 
0    TO    854 


LDtLPD,JJ),K,LDENtLX,K),JJJ) 

EQ.PARAMt 1).AND.LDEN(LX,2) . EQ.P ARAM t 2 ) . A 
EO. PARAMO)  )     GO    TO    840 
LDtLPD,JJ) ,4,6,4) 


LDtl, JJ) ,4,7,4) 
855,60 


1 

(KK) 

I,LK1, JS,  t  tLNUM(K,Kl),Kl=l,3),K=l,KK),IS 


(KK) 
H ( KKK ) 

I,LK1,JS, ( (LNUMtK,Kl) ,K1=1,3) ,K=1,KK) , t  t 
1 ,3) ,K=1 ,KKK), ISLOP 
W)232,901,232 

9  02,904 

DP 


,932 
OP 


,932 


T,4,PARTLNtLPN,K),4) 

,60,961 

LNtLPN,K) ,4,4,4) 

LNtLPN,K) ,4,9,4) 


T,4, PARTLDtLPD,K),4) 

,60,963 

LDtLPD,K)  ,4,4,4) 

LDtLPD,K) ,4,9,4) 
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931  CONTINUE 

932  IF(N9.E0.N7)     30    TO    909 
GO    TO    24 

909       IF(WHERE-1 )24,930,930 
930       CALL     SENS 

24       CONTINUE 
10  24    CONTINUE 

IPRIN-IPRIN-1 

IF( IPRIN)13,13,89    " 

5  LIC=LIC+1 
WRITE    (6,163) 
DO    12    N1=1,M9 
IF(LIC-IABS( ILOOP(Nl) ) )4,4,12 

4  DO    11    N7=1,MVI 

IF(M0DNM(N1)-I ABS(NAME(N7,1) ) 111,6,11 

6  IF( ILOCP(Nl) ) 8,60,7 

7  ZZ(N7)=VALM0D(N1,LIC) 
GOTO    9 

8  ZZ(N7)=VALM0D( Nl , 1 ) +FLOAT ( L IC-1 )* VALMODt Nl ,3) 

9  NM=NAME(N7,  D/100 

WRITE  (6,129)  (NAM(N7,K) ,K=1,3) ,ZZ(N7) 
IF( (NM+3)*(NM-4)*(NM-5) ) 12, 10, 12 

10  ZZ(N7)=1./ZZ(M7) 
GOTO  12 

11  CONTINUE 
WRITE  (6,130) 
IERR=1 
MVI1=MVI+1 

DO  15  N7=MVI1,  M 

IF( M0DNM(N1)-I ABS(NAME(N7,1) )) 15, 14,15 

14  WRITE  (6,133)  MODNM(Nl) 
GOTO  12 

15  CONTINUE 

WRITE  (6,132)  MODNM(Nl) 

12  CONTINUE 

IF{ IERR)60,13,  66 

13  IF( ILOOPfl  )-1300)17,18,60 
18    WRITE  (6,101) 

17    DO  400  Nl=l,30 
SDEN(N1 )=0. 
DO  400  K2=1,MTRN1 

400  SNUM(N1,K2)=0. 

DO  603  K2=1,MTRN1 
603   SNUM( 16,K2)=N^BR(K2) 

SDEN( 16)=1. 

1=0 

DO  1412  K2=1,MDD 

N8  =  l 

N7=LIM(K2) 

IF(N7)60, 1412, 397 
397       DO   412    N9=1,N7 

1  =  1  +  1 

ISLOP=0 

VL00P=1. 

LK1  =  L00P(I  ) 

LK2=IABS(LK1  ) 

SGN=LK1/LK2 

401  N1  =  N8 

402  IF(LK2-ITW0(N1 ) )    403,405,405 

403  IF(N1-N8)60,404,406 

404  N8=N8+1 
GOTO    401 

405  ISLOP=I SL0P+IST(N1) 
VL0CP=VL00P*ZZ(N1) 
LK2=LK2-ITW0(N1 ) 

406  N1=N1+1 

IF (Nl-M VI- 1)402, 407, 60 

407  JS=16-ISL0P 
IF(MDD-K2)     60,410,408 

40  8   SNUM(JS,K2)=SNUM( J S , K2 ) +VLOOP*SGN 

GOTO  412 
410   SDEN( JS)=SDEN( J S )+ VLOOP*SGN 
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412       CONTINUE 
1412    CONTINUE 

DO    69    K=1,MTRN1 

JZ(K)=1 

AZ(K)=0. 

DO    68    J  =  l, 30 

IF{SNUM(J,K))67,68t67 
67  IF(AZ(K) )201 ,202 ,231 

202       AZ(K)=-SNUM( J,K) 
201     JZ(K)=J 

68  CONTINUE 

69  CONTINUE 
SIGNC=0. 
JI=0 

DO    71    J=l,30 
IF(SDEN{ J) )97, 71,97 

97  IF(SIGNC) 70,98,70 

98  SIGNC=SDEN(J ) 
73  JI=J 

71  CONTINUE 
J  =  0 

DO   73  K=1,MTRN1 
IFUZ(K)-J)  73,73,72 

72  J=JZ(K) 

73  CONTINUE 
IF(JI-J)   75,75,74 

74  J=JI 

75    DO  83  K=1,MTRN1 

XG(K)=0. 

IF(SIGNC)90,91 ,90 
90    XG(K)=AZ(K) /SIGNC 

91  WRITE  (6,136)  XG(K) 

WRITE  (6,135)  ( ITAG(K, I) ,1=1,8) 

NU  =  0 

IO(K)=0 

DO   82  L=1,J 

ANUM=NU 

IF(ABS(SNUM(L,K) )*(ANUM-.5) ) 8  0, 82, 81 

80  IO(K)=J-L 

81  NU=NU+1 
IP(NU,K)=J-L 

SNUM(NU,K) =-SNUM(L,K)  /AZ(K) 

WRITE     (6,137)     SNUMINU,K) ,IP(NU,K) 

82  CONTINUE 
NUM(K)=NU 
IF(NUM(K))     60,92,83 

92  NUM(K)=1 

WRITE     (6,137)        SNUM(l.K) ,IP(1,K) 
WRITE     (6,141) 

83  CONTINUE 
WRITE     (6,138) 
NUMZ=0 

ID=0 

DO       86    K=1,J 

ANUM=NUMZ 

IF(ABS(SDEN(K)  )*(ANUM-,5) ) 84, 86,  8 5 

84  ID=J-K 

85  NUMZ=NUMZ+1 
IOD(NUMZ)=J-K 
SDEN(NUMZ)=SDEN(K) /SIGNC 

WRITE     (6,137)     SDEN(NUMZ),IOD(NUMZ) 

86  CONTINUE 
IF(NUMZ)    60,9^,87 

94    NUMZ=1 

WRITE     (6,137)        SDEN(l) , I OD ( 1 ) 
WRITE    (6,141) 
87  IFILI-LIO150,  150,  152 

153       LIC=3 

152  IF(MTRN1-MTRN2 )96,96,95 

95  IFCMTRN1-MTRN3 ) 88,88,  153 

153  IF(LIC) 60,89,5 

96  IND=5 
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GO  TO  900 

88  IND  =  6 

GO  TO  900 
60    WRITE  (6,134) 
66    LIC=0 

89  IND=1 
9  00  P c  TURN 

100  F0RMAT(//1X,16HIDENTICALLY    ZERO) 

101  F0RMAT(///1X,16HN0MINAL    SOLUTION) 

129  FORMAT! /1X  ,  1  OH  MODI F YI NG     ,A1,2I1,4HT0    ,E11.4) 

130  FORMATt //1X, 26HERR0R    IN    ATTEMPT    TO    MODIFY) 

132  FORMAT! IX, 12HBRANCH    CODED, 15, 22H — NOT    IN    ORIGINAL    NETW 
■j      3H1 P  K  ) 

133  FORMAT! IX, 12HBRANCH  C0DED,I5,22H — CANNOT  MODIFY  INDEPE 
1,13HNDENT    SOURCES) 

134  FORMAT! //1X, 24HMACH INE    LANGUAGE    FAILURE! 

135  F0RMAT!/1X,24HNUMERAT0R  POLYNOMIAL  OF  , 2A 1 , 21 1 , 1H/, 2A1 
1,211) 

136  FORMAT! //1X,4HK    =     ,E11.4) 

137    FORMATQX,  E11.4,  10H   TIMES    S**,I3) 
138       FORMAT! //1X , 22 HDE N3MI NATOR    POLYNOMIAL) 

141    FORMAT!/, IX, 31HCHECK    YOUR    CIRCUIT    CODING    LIST., 
+/,1X,26HTHIS    ANSWER    IS    DEGENERATE.) 
160       F0RMAT(/3X,28HL0C         CODE  LOOP    PRODUCT) 

162  F0RMATU7X,  Al,  13H(     1     )     S    **       0) 

163  FORMAT!/) 

164  FORMAT!'  ',» WARNING:  NUMBER  OF  UNCANCELLED  LOOP  PRODU* 
1,«CTS  EXCEEDS  STORAGE  AVAILABLE  (200).  «,/«  ','YOU  MU« 
2,' ST  CHANGE  DIMENSION  STATEMENTS  IN  NAS4  AND  SENS  SUB* 
3, 'ROUTINES  FOR  VARIABLES  -  ',//«  «  , ' PARTLN  ,  EXPNS ,  NFN« 
4, ■ , PART LD, EXPDS, SIGD,NF D, SIGN, RNUMP,RNUMS,RDENS,RDENP« 
5) 

STOP 
END 

PACKING  SUBROUTINE 

SUBROUTINE  PAC K! A , I ,B , J) 

L0GICAL*1  A! 1) ,B!1) 

A(I)=B! J) 

RETURN 

END 


SUBROUTINE  SENS 

COMPLEX* 8  CFRE0,CTEMP,CTEMP1,CTEMPB 

INTEGER  EXPNS, S IGN , EXPDS , S IGD, WHERE , PAR AM , PARTLD, 
1 PARTLN 

DIMENSION  RNUMP! 2  0  0)  ,RNUMS!2  00)  ,RDENS(2  00)  ,RDENP!200), 
1  ITEMP!23), IEXPN! 50) , IEXPD! 50) , PARTI! 50)  ,PART2(50) 

C0MM0N/SEN/PARAM!3) , SENSFL , FRM AX , FRMIN 

COMMON  IND,NASA,MTRN1,MTRN2,MTRN3,MD,IJTAG!15,2) , I 
1PRIN,Z!42) ,LL, M0DNM(23) ,VALMOD( 23, 10), I  LOOP (23) ,LI,LIC 
2,M9,ZZ!23)  , NAM! 23,3  )  , I  TAG! 5 ,8)  ,NMBR(6), 
3M,SNUM( 3  0, 5) ,SDEN( 30) , I0D(30) ,10(5) ,NUM(5) 
5,XG(5),IP(30,5), ID,NUMZ,DMY(23) 

4, NAME! 23,2),MVI,IST!40),LIM(6) , I  TWO (23) , LOOP (10000) 
6,WHERE,PARTLN( 200, 23 ) , E XPN S 
7(200),NFN(200),PARTLD(200,23), EXPDS (200), 
8SIGD(200),NFD( 200)  ,LPN  ,  LPD , SI GN( 200 ) 

DATA  INEG/1-    '/,IBLNK/'     V , PART  1/50*0. 0/, PART2/50 
1*0.0/ 

DO  2  K=1,MVI 

ITEMP(K)=NAM(K,1 ) 

DO  2  L=2,4 

CALL    PACKUTEMP(K)  ,L,0,4) 

CONT  INUE 

IF(WHERE-1)60,200,19 
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19 


200 


205 


210 

220 
225 
243 


25D 

280 
285 

290 

300 

310 
315 
316 

317 

320 
400 


402 


405 


410 

420 
425 

440 


450 


LPN 


KP 


CALL  PACK(PARTLD 

CALL  PACKtPARTLD 

NFD( 1)=1 

GO  TO  400 

DO  320  I =1 

RNUMPt I } =1 

RNUMSt I )=1 

KP=NFN(I ) 

DO  320  J=l 

ICOMP1=0 

ICOMP2=0 

ICOMP3=0 

IC0MP4=0 

CALL  PACKUC0MP1 

CALL  PACK(IC0^P2 

CALL    PACK(IC0*iP3 

KL  =  0 

KL=KL+1 

IF( IC0MP1.EQ.ITE 

GO  TO  225 

IF{  IC0MP2.EQ.NAM 

GO  TO  225 

IF( IC0MP3.E0.NAM 

IF(KL.LT.MVI)  GO 

CALL  PACKUC0MP4 

IC0MP4=IC0MP4-H 

GO  TO  (250,60, 60 

RNUMPt I )=RNUMP( I 

GO  TO  285 

RNUMP( I )=0 

RNUMSt I )=RNUMSt I 

GO  TO  316 

RNUMPt I )=RNUMP (I 

GO  TO  315 

RNUMP(  I  )=RNUMP (  I 

RNUMPt I ) =-RNUMP( 

GO  TO  315 

RNUMPC I )=0 

RNUMSt I ) =RNUMS (I 

IFtSIGNt IJ.EO. IN 

GO  TO  320 

RNUMPt I ) =-RNUMP( 

RNUMSt I )=-RNU^S( 

CONTINUE 

GO  TO  30 

DO  2520  1=1, LPD 

KP  =  NFD( I  ) 

RDENPt I )=1 

RDENSt I )=1 

IFtl.NE.l) 

RDENPt I )=0 

GO  TO  2520 

DO  520  J=l 

IC0MP1=0 

IC0MP2=0 

IC0MP3=0 

IC0MP4=0 

CALL  PACKt IC0MP1 

CALL  PACKUC0MP2 

CALL  PACKt IC0^P3 

KL=0 

KL=KL+1 

IF  (  ICCWP1.EQ.ITE 

GO    TO    425 

IF (  IC0MP2.E0.NAM 

GO  TO  425 

IFI IC0MP3.EQ.NAM 

IFtKL.LT.MVI  )  GO 

CALL  PACKt IC0^P4 

IC0MP4=IC0MP4f 1 

GO    TO    (450,60,60 

RDENPt I )=RDENP(  I 


(1,1)  ,3,1,4) 
tl, 1) ,4,4,4) 


,1 ,PARTLN{ I,J),1) 
,4,PARTLNt I , J) ,2) 
,4,PARTLNt I, J), 3) 


MP(KL) )    GO    TO    210 

tKL,2) )     GO    TO    220 

tKL,3) )     GO    TO    240 

TO    205 
,4,PARTLNt I, J) ,4) 


,2  85,280,60,290,300,60, 310), ICOMP 4 
)/ZtKL) 


)/ZtKL) 
)*Z(KL) 


)*Z(KL )**2 
I) 


)*Z(KL) 

EG)    GO    TO    317 


I) 
I  ) 


GO    TO    402 


KP 


,l,PARTLDt I, J), 1) 
,4,PARTLD(I ,J) ,2) 
,4,PARTLD( I, J) ,3) 


MP(KL) )    GO    TO    410 

tKL,2) )     GO    TO    420 

tKL,3) )     GO    TO   440 

TO    405 
,4,PARTLDt I, J) ,4) 


,485,480,60,490,500,60,510) ,IC0MP4 
J/ZIKL) 
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GO  TO  485 
480   RDENP(I)=0.0 
485   RDENSII  )=RDENS  (I  )/Z(KL) 

GO  TO  516 
490   RDENPl I)=RDENP(I )*Z(KL) 

GO  TO  515 
500   RDENP( I  )=RDEND ( I  )*Z(KL)**2 

RDENPt I )=-RDEMP( I) 

GO  TO  515 
510   RDENP(I)=0.0 

515  RDENSU  )=RDENS(  I)*Z(KL) 

516  IF(SIGD( I) .EO. INEG)  GO  TO  517 
GO  TO  520 

517  RDENPfl )=-RDENP(I) 
RDENSC I ) =-RDENS(I ) 

520   CONTINUE 
2520  CONTINUE 

J  =  0 

DO  610  K=1,LPD 

DO  610  1=1 tLPN 

J  =  J  +  1 

IEXPNt J)=EXPNS ( I )+EXPDS(K) 

IF( ABS(RDENS(K)*RNUMP( I ) -RNUMS ( I J*RDENP (K)  ).LE.1.0E-75 
1)  GO  TO  1520 

GO  TO  1521 

1520  J=J-1 

GO  TO  610 

1521  IF(J.EO.l)  GO  TO  1620 
JMIN1=J-1 

DO  1610  KL=1,JMIN1 

IF(IEXPN(J).EQ.IEXPN(KL) )  GO  TO  608 

GO    TO    1610 
608       PARTI  (KL)=PART1(KU+RDENS(K)*RNUMP(  I  J-RNUMSC  I  )*RDENP(K 
1) 

J  =  J-1 

GO    TO    610 
1610    CONTINUE 

16  20    PARTI (J )=RDENS(K)*RNUMP(I)-RNUMS(I)*RDENP(K) 
610       CONTINUE 

IF(J.LE.l)     GO    TO    640 

JJ  =  1 

631  K=JJ+1 

IF(A8S(PART1(J J) ). LE.1.0E-75)     GO    TO    636 

632  IF( lEXPN(JJ) .GT.IEXPN(K) )    GO    TO    634 
TEMP1=PART1 (JJ ) 

TEMP2=IEXPN(  JJ  ) 

PART1(JJ)=PART1(K) 

IEXPN( JJ)=IEXPN{ K) 

PART1(K)=TEMP1 

IEXPN(K)=TEMP2 
634       K=K+1 

IF(K.GT.J)    GO    TO    636 

GO    TO    632 
636      IF( JJ.EO.C J-l) J    GO    TO    640 

JJ=JJ+1 
638       GO    TO    631 

640  WRITE(6,900)     (  I  TAG < It K)  ,K=1 ,8)  » (PARAMC K)  ,K=1 ,3 ) 
IF(J)60,641,642 

641  WRITE(6,906) 
GO    TO    643 

642  SET=0 

DO    710    K=1,J 

IFUBS(PARTHK)  )  .LT.1.0E-75)    GO    TO    710 
SET=1 
700    WRITE(6,901)     PARTI ( K) , IEXPNl K) 
710       CONTINUE 

IF(SET.EO.O)    dRITE(6,906) 

643  J1=J 
J  =  0 

DO  611  K=1,LPD 
DO  611  1=1 ,LPD 
J  =  J  +  1 
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IEXPD( J)=EXPDS ( I )+EXPDS(K) 

IF(ABS(RDENS(K)*RDENS( I)) .LE.l.OE-75)    GO    TO    1711 
GO    TO    1712 

1711  J  =  J-1 

GO    TO    611 

1712  IF(J.EO.l)    GO    TO    1612 
JMIN1=J-1 

DO    1611     KL=ltJMINl 

IF(IEXPDU).EQ. IEXPD(KL) )     GO    TO    609 
GO    TO    1611 
609       PART2(KL)=PART2(KL)+RDENS(K)*RDENS(I ) 
J  =  J-1 
GO    TO    611 

1611  CONTINUE 

1612  PART2C J)=RDENS(K)*RDENS< I) 
611       CONTINUE 

IF(  J-D60.749,  740 

740  JJ=1 

741  K=JJ+1 

IF(ABS(PART2(J J)  ) . LE. 1 . OE-75 J     GO    TO    746 

742  IF( IEXPD(JJ).GT.IEXPDIK)  )     GO    TO    744 
TEMP1=PART2( JJ ) 

TFMP2=IEXPD( JJ  ) 

PART2( JJ)=PART2(K) 

IEXPDt JJ)=IEXPD(K) 

PART2(K)=TEMP1 

IEXPD{K) =TEMP2 
744       K=K+1 

IF(K.GT.J)    GO    TO    746 

GO  TO  742 
746   IF(JJ.EQ.(J-1)  )  GO  TO  749 

JJ=JJ+1 

748  GO  TO  741 

749  WRITE(6t902) 
DO  750  K=l ,  J 

IF(ABS(PART2(K) ) .LT.l. OE-75)  GO  TO  750 
WRITE(6,901)  PART2(K),  IEXPD(K) 

750  CONTINUE 

IF(FRMAX.EQ.O. 0. AND . FRMI N. EQ . 0. 0)    GO    TO    30 

FRINC=(FRMAX-FRMIN  )/ 199.0 

FREO=FRMIN 

WRITE(6,904) 

DO    763     KK=1,200 

FREO 1=2. 0*3. 14 1592 6535* F RE Q 

CFRE0=CMPLX(0.0»FREQ1) 

CTEMP=(0.0t0.3  ) 

DO    755    K=1,J1 

CTEMPl=CMPLX(PARTl(K)t 0.0) 

IFUEXPN(K)  .LT.O)     GO    TO    754 

CTEMP=CTEMP+CTEMP1*CFREQ**IEXPN(K) 

GO    TO    755 
7  54    CTEMP=CTEMP+CTEMP1/CFREQ**IIABS(IEXPN(K) ) ) 
755      CONTINUE 

CTEMPB={0. 0,0.0  ) 

DO    760    K=ltJ 

CTEMP1=CMPLX(PART2(K) ,0.0) 

IFUEXPD(K)  .LT.O)    GO   TO    759 

CTEMPB=CTEMPB+CTEMP1*CFREQ**IEXPD(K)   . 

GO    TO    760 
7  59    CTEMPB=CTEMPB+CTEMP1/CFREQ**(I ABS(  IEXPD(K  )  )) 
760      CONTINUE 

IF(CABS(CTEMPB ) . LT . 1. OE-74 )    GO    TO    1762 

CTEMP=CTEMP/CTEMPB 

RMAGN=CABS(CTEMP) 

ANGLE=ATAN2(AIMAG(CTEMP) , RE AL ( CTEMP ) ) 

762  WRITE(6,905)     FRE Q, RMAG N, ANGLE 
1762     FREO=FREO+FRINC 

763  CONTINUE 
GO    TO    30 

900       FORMATt'     • , / IX , 3 1HNUMERAT0R    POLYNOMIAL     :     SENSITIV 
1  ,'ITY    OF    «  ,2A1,2I1 , 
21H/,2Al,2Ut2X,4HT0       ,A1,2U) 
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901  FORMAT («  • ,  / , IX » El  1 .4,  10H  TIMES  S**tI31 

902  FORMAT (»  » t /// / , IX , 27HDEN0M IN ATOR  POLYNOMIAL 
1, 'ENSITIVITY'  ) 

903  FORMAT (•  ■,//,•••,■ MACHINE  FAILURE') 

904  FORMATC  »,////•  ■  ,  ■  EV  ALUAT  I  ON  OF  SENSITIVITY 
1///,'  ',  'FREQUENCY' 

2  aOX,' MAGNITUDE'  , 10X, «PHASE(  IN  RADIANS)  '  t /,  ' 

905  FORMATC  •  ,E  11  .  4  ,  8X  ,E1  1  .4  ,  8X  ,  El  1  .  4  ) 

906  FORMAT (■  '»//♦'  ' ,  ' I  DENT  ICALLY  ZERO',//,'  f) 
60   WRITE{6,903) 

30    RETURN 
END 
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